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Abstract 



^ , Level shift operators describe the second order displacement of eigenvalues under per- 

• turbation. They play a central role in resonance theory and ergodic theory of open quantum 

systems at positive temperatures. 

We exhibit intrinsic properties of level shift operators, properties which stem from the 
structure of open quantum systems at positive temperatures and which are common to 
all such systems. They determine the geometry of resonances bifurcating from eigenvalues 
of positive temperature Hamiltonians and they relate the Gibbs state, the kernel of level 
shift operators, and zero energy resonances. 
, We show that degeneracy of energy levels of the small part of the open quantum system 

O |i causes the Fermi Golden Rule Condition to be violated and we analyze ergodic properties 

1 j of such systems. 

C3 



1 Introduction and main results 



Level shift operators emerge naturally in the context of perturbation theory of (embedded) 
eigenvalues, where they govern the shifts of levels (resonances) at second order in perturbation. 
They play a central role in many recent works on ergodic properties of open quantum systems 
at positive temperature [20, 21, 9, 14, 26, 17, 18, 19, 23, 28, 2, 1]. The dynamics of such systems 
is an automorphism group of the algebra of observables generated by an operator L\ = Lq + XI, 
where the selfadjoint Lq describes the free dynamics of two (or more) uncoupled subsystems, 
A £ K is a coupling constant, and I is an interaction operator. 

Ergodic properties are encoded in the spectrum of L\. If the system has an equilibrium 
state at positive temperature 1/(3 one shows that L\Qp t \ = and that if Ker L\ = CQp^x 
then any state initially close to equilibrium approaches the equilibrium state in the limit of 
large times. (We do not address the question of mode or speed of the return to equilibrium in 
this outline). 

It follows from the algebraic structure of quantum systems at positive temperatures that the 
operator Lq has necessarily a degenerate kernel whose elements are in one-to-one correspondence 
with invariant states of the uncoupled system. In order to prove return to equilibrium one needs 
to show that the degeneracy of the eigenvalue zero, which is embedded in continuous spectrum, 
is lifted under perturbation: dim Ker La = 1 for A ^ 0. This has been proven for several 

'Supported by a CRM-ISM postdoctoral fellowship and by McGill University; merkli@math.mcgill.ca; 
http: / /www. math. mcgill.ca/~merkli/ 
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concrete models [20, 21, 9, 26, 14, 18, 2]. The fate of embedded eigenvalues under perturbation 
can be described by spectral resonance theory if the system has certain deformation analyticity 
properties, [20, 21, 9, 2], or, for less regular systems, by a Mourre theory [14] or by a positive 
commutator theory, [26, 18]. 

A core strategy common to these methods is to reduce the spectral analysis of the operator 
L\ around the origin to that of a reduced operator acting on a smaller Hilbert space (which 
is finite-dimensional in all works cited above). This procedure is sometimes called 'integrating 
out degrees of freedom'. For deformation analytic systems it can be implemented by applying 
the so-called Feshbach map F [8] to a suitably deformed operator L\(r), where t e C is the 
deformation parameter 1 . The Feshbach map has an isospectrality property, implying that the 
kernels of F(L\(t)) and of L\ are isomorphic. An expansion in the coupling constant A gives 

F(L A (r)) = -A 2 A + O(A 3 ), (1) 

where the operator Ao is independent of the deformation parameter r. The property of return 
to equilibrium follows if A has simple kernel because then (1) and the isospectrality of the 
Feshbach map imply that for small A ^ 0, dimKerL^ < 1, and since Lx^lp^x = one must 
have KerL\ = CQp^x. The operator A is called the level shift operator (associated to the 
eigenvalue zero of La)- For a more detailed description of the emergence of level shift operators 
in perturbation theory we refer to the works cited above, and also to [12]. 

Level shift operators are equally important in the study of systems far from equilibrium, 
where the system does not have an equilibrium state, for instance when several thermal reser- 
voirs at different temperatures are coupled, [23, 28, 1], or when the small system does not 
admit an equilibrium state, [17, 19]. In what follows we discuss the former case. The role of 
the equilibrium state is now played by a reference state tpx, e.g. the product state of the small 
system and the reservoirs in equilibria at different temperatures. An interaction operator W 
can be chosen such that the Heisenberg dynamics of the system is generated by the operator 
Kx = L + XW, satisfying Kxipx — 0. Unlike / in the situation of systems close to equilibrium, 
the operator W cannot be chosen to be normal. A detailed spectral analysis of operators of this 
type (called 'C-Liouville operators') is carried out in [28]. One obtains the level shift operator 
A associated to the eigenvalue zero of Kx by an expansion of F(Kx(t)) in A, as in the situation 
above. In the context of systems far from equilibrium, a dynamical resonance theory shows 
that if A has simple kernel then the system possesses a unique time-asymptotic limit state (for 
small A), which is a non-equilibrium stationary state. 

We should also point out that level shift operators have a dynamical interpretation as 
"Davies generators" of the reduced dynamics in the van Hove limit, [14, 13]. 

In the present paper we examine properties of level shift operators which do not depend on 
particularities of the system in question, but which originate from the structure common to all 
open quantum systems at positive temperatures. One of the main ingredients determining this 
structure is the Tomita-Takesaki theory of von Neumann algebras. 

We describe the geometry of resonances bifurcating from real eigenvalues of Liouville oper- 
ators in Theorem 1.1. Part (c) of that theorem examines the role of the degeneracy of energy 
levels of the small system. The interplay between the Gibbs state of the small system and the 
kernel of the level shift operator is described in Theorem 1.2. Some of these intrinsic proper- 
ties we exhibit here have been observed in the analysis of specific systems carried out in the 
references given above. 

Apart from an analysis of the structure of level shift operators we study in Section 3 the 
dynamics of systems where a Bosonic heat reservoir is coupled to a small system whose Hamil- 
tonian has degenerate eigenvalues. For such systems the so-called Fermi Golden Rule Condition 

lr rhc operator L\ itself is typically not in the domain of the map F but L\(t) is for r K. Arguments 
similar to the ones we give here also work for systems which are not deformation analytic, but they need a 
technically more elaborate presentation. 
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is violated due to the fact that the interaction "cannot couple the Bohr frequency zero" (arising 
from the degenerate levels) to the zero energy reservoir modes in an effective way. We prove 
return to equilibrium for such systems by taking into account higher order corrections in the 
perturbative spectral analysis of the operator L\ . A consequence of the degeneracy of energy 
levels of the small system is that the approach to equilibrium is still exponentially fast but with 
relaxation time of 0(A -4 ) as opposed to the shorter relaxation time 0(A -2 ) for systems with 
non-degenerate spectrum. 

The organization of this paper is as follows. In Section 1.1 we summarize some facts about 
open quantum systems. Our main results concerning the structure of level shift operators are 
Theorems 1.1, 1.2 and 1.4, given in Section 1.2. In Section 2 we present examples of concrete 
models to which our results apply. We analyze systems with Hamiltonians having degenerate 
eigenvalues in Section 3. Section 4 contains proofs. 

1.1 Open quantum systems 

A detailed description of open quantum systems can be found in [10], in the above-mentioned 
works, and also in [22, 27]. 

Consider a quantum system possessing finitely many degrees of freedom, like a single particle 
or a molecule, or a system with finitely many energy levels. We denote by JC the Hilbcrt space 
of pure states of this "small system" and by H its Hamiltonian. We allow the case dim/C = oo 
but require that the Gibbs state ^ p exists, i.e., that e _/3H is trace class on JC, for some inverse 
temperature < (3 < oo. In the above example the particle or molecule must thus be confined, 
e.g. by a potential. 

We view the Gibbs state ^asa vector in the Hilbert space Hs — JC ® JC. The algebra of 
all bounded operators, B(JC), contains the observables of the small system. It is represented 
on Hs as the von Neumann algebra Wis = B(JC) <g> i/c C B(Hs), and the dynamics on Wis is 
implemented by 1 i— > e ltis ^4e~ ltis , where 



is called the standard Liouville operator of the small system. The Gibbs vector typ is cyclic 
and separating for Wis- We denote the modular conjugation associated to the pair (Wls,^p) by 
Js- The standard Liouville operator Ls satisfies the relations Lg^fp — and JsLsJs = —Ls- 

In models of systems close to equilibrium, the small system is placed in a (single) envi- 
ronment (reservoir) modeled by a "large" quantum system having infinitely many degrees of 
freedom. A common example is a spatially infinitely extended ideal quantum gas (of Bosons or 
Fermions). We assume that the reservoir has an equilibrium state (for some inverse temperature 
< [3 < oo) which is represented by the vector $>p in the reservoir Hilbert space Kr. 

Observables of the reservoir are operators belonging to (or affiliated with) a von Neumann 
algebra 9JIr c £>(7Yr). Their dynamics is given by a group of automorphisms of Win, t i— > 
grf-t-R Ae~ ltLn , generated by a selfadjoint standard Liouville operator Lr. Being a KMS vector 
w.r.t. this dynamics, <f>p is cyclic and separating for $Wr. We denote the modular conjugation 
associated to (Wln,$p) by Jr. The operator Lr has the properties JrLrJr = — Lr and 
= 0. 

The von Neumann algebra Wl = Wis ® 9#R, acting on the Hilbert space H — Hs <8> Hr, 
contains the observables of the combined system. Elements in this algebra evolve according to 
the group of automorphisms of Wl generated by the selfadjoint operator 



L s = H ® 1 K - I*; ® H 



(2) 



Lo — Ls + Lr. 



(3) 



The vector 



(4) 
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defines the equilibrium state w.r.t. this dynamics, at inverse temperature (3. 

To describe the interacting dynamics we introduce the map it which sends linear operators 
on K ® Hr to linear operators on K. <E> K. <£> Hr according to 

tt(A <g> B) = A <g> t K <g> S. (5) 

The interaction between the small system and the reservoir is specified by a selfadjoint operator 
v on /C <g> such that 

V = 7r0) (6) 
is a selfadjoint operator affiliated with 9JI. Let J = Js <£> Jr. We assume that 

Bl L + XV is essentially selfadjoint on T>(L ) n ©(F) and i + XV — XJVJ is essentially 
selfadjoint on T>(L ) n 2?(V) D £>(JVJ), where A is a real coupling constant. 

If assumption Bl holds then the selfadjoint operator 

L x = L + XI, (7) 

where 

I = V- JVJ, (8) 

generates a group of automorphisms a\ — e ltLx ■ e~ ltLx of 9JI (see e.g. [16], Theorem 3.5). We 
are interested in interactions for which the coupled dynamics a\ admits an equilibrium state. 
It is known that the condition 

B2 ty, i0 € P( C -/3(io+AV)/2) 

implies that 

filX ~ \\c-^ L o+^V)/2 n/3 Q \\ ^> 

is a (/3, a^)-KMS state, and that the following properties hold: JL\J = —L\, L\Q,p x — 0, and 
lim^^o 11^/3, a — ^/3,oll = (sec e.g. [16], Theorem 5.5). 

In models for systems far from equilibrium the small system is coupled to several reservoirs. 
The Hilbert space of the small system plus R > 2 reservoirs is TC = Hs ® Hr ® ■ ■ ■ <8> Hr and the 
non-interacting dynamics of the algebra 971 = 50ts ® 9^R ® • • • 5Hr is generated by the selfadjoint 
standard free Liouville operator 

R 

L = L s + Y,Ln,r, (10) 

r=l 

where iR !r acts non-trivially only on the r-th reservoir space. The interaction is determined 
by an operator 

R 

K = 5>( Vr ), (ii) 

r=l 

where w r is the selfadjoint operator on K ® Hr representing the coupling between the small 
system and reservoir r. We understand that n(v r ) acts trivially on all reservoir Hilbert spaces 
except the r-th one. As mentioned in the introduction, the role of the equilibrium state is now 
played by a reference state ipx € H. We assume the following. 

C The interacting dynamics of 9Jt is generated by the operator Kx = L n + XW, where W is 
a bounded, linear (generally non-symmetric) operator on H, s.t. Kxipx — for all A in a 
neighbourhood of zero. 

Remark. For systems with bosonic heat reservoirs the operator W is not bounded. Such 
systems are considered in [28]. The development of a general theory for unbounded non- 
symmetric W is a technically intricate affair, we restrict our attention in this note to systems 
satisfying condition C (although we give a more general result in Theorem 4.2). 
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1.2 Main results 



1.2.1 Systems close to equilibrium 

We assume that the basic assumptions Bl and B2 are satisfied. The spectral projection onto an 
eigenvalue e of L$ is denoted by XL s =e and Pr = I^X^I denotes the orthogonal projection 
onto C$0. Set P e — XL s =e <8> Pr and P e = 11 — P e . We allow the case dimP e = oo. Denote 
by L x the restriction of L\ to RanP e , i.e., L x := P e L\P e [ Ran p e - Consider the following 
condition: 

Al e P e I and IP e are bounded operators. 

If condition Al e holds then we define the family of bounded operators 



for real e ^ 0. Note that for e > we have ImA e (e) = eP e IP e [(L e Q - e) 2 + e 2 ] _1 P e /P e > 0, so 
the numerical range, and hence the spectrum of A e (e), lie in the closed lower complex plane. If 
(12) has a limit as e J. (in the weak sense on a dense domain) then we call this limit the level 
shift operator associated to the eigenvalue e, and write it as 



The projection Pr has rank one so it is natural to identify RanP e with R&nxL s =e- In this 
sense we view the operators (12) and (13) as operators acting on Ranxi s =e C JC <g> JC. 

Theorem 1.1 

(a) Let e be an eigenvalue of L$. Al e holds if and only i/Al_ e holds. If the assumptions 
Al± e hold then JsA e (e)Js = — A_ e (e). 

(b) Assume AIq. Suppose the spectrum of the Hamiltonian H consists of simple eigenvalues 
and denote the eigenvectors by tpi. Then Ao(e) = iro(e), where To(e) is a self adjoint 
positive definite operator on RanP which has real matrix elements in the basis {ipi <g> ipi}. 
If A exists then the same statements are true for that operator. 

(c) If H has degenerate eigenvalues then neither the real nor the imaginary part of Ao(e) 
vanish, in general. The matrix elements of A (e) in the basis {ipi ® ipj} are not purely 
real nor are they purely imaginary, in general. If A exists then the same statements are 
true for that operator. 

Part (a) of the theorem shows that <r(A e (e)) = — (j(A_ e (e)). In particular, the spectrum of 
Ao(e) is invariant under reflection at the imaginary axis. Part (b) shows that <r(Ao(e)) lies on 
the negative imaginary axis if H has simple spectrum. Part (c) says that if H has degenerate 
spectrum then cr(Ao(e)) can have nonzero real part. In the context of a translation-analytic 
model of an iV- level system coupled to a bosonic (or fermionic) heat reservoir [20] (see also [28]) 
Theorem 1.1 shows that the set of all resonances is invariant under reflection at the imaginary 
axis, and that resonances bifurcating from the origin stay on the imaginary axis while wandering 
into the lower complex plane if H is non-degenerate (see also (1), (61) and (62)). 

We present a proof of assertions (a) and (b) of Theorem 1.1 in Section 4.1 below. In Section 
3 we give examples illustrating statement (c). 

Our next result concerns the interplay between the KMS state f2/3,A, (9), and the level shift 
operator for e = 0, A . We introduce the following assumptions. 

A2 A i ► P IPo(L\ — ie) _1 Po is continuous at A = as a map from R to the bounded 
operators on H, for nonzero (small) e. (We write here L\ instead of L x .) 



A e (e) = P e IP e (L e - e - ie)- 1 P e /P c 



(12) 



A e = P e IP e (L e - e - iO+)- 1 P e /P e . 



(13) 
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A3 A i ► is diffcrentiable at A = as a map from R to H. 
Theorem 1.2 Assume that Conditions Al , A2 and A3 hold. Then we have 

lim PoIPo(L - ie^PoIPon^o = P IP d A | A=0 %A- (14) 
In particular, if Ao exists and if PqIPo — then Ao^p = 0. 

Even though wc do not assume that A e (e) converges, equation (14) shows that A e (e)fi / 3 i o 
has a limit as e — ► 0, regardless of the sign of e. By subtracting from (14) that same equation 
with e replaced by — e we get 

P /P S(L ) fV^otyj.o = 0, (15) 

and by adding the two equations we obtain 

P IP P.V.(Lo)- 1 Po/PoV^o = PoIPo 9a|a=o%a, (16) 

where 8{x) = lim e |o \ x ^+^ * s ^ ne Dirac delta distribution and P.V..T -1 = lim e |o x -i+ e -i is the 
principal value distribution. (The limits are understood in the strong sense.) 

Remarks. 1) If A exists then (15) shows that the Gibbs state ^ p of the small system 
belongs to the kernel of Im Ao, c.f. (4). If PqIPo = 0, then (14) implies the second statement 
of the theorem. 

2) If dim/C < oo, PoIPo = 0, and if A exists, a characterization of Ker(ImA ) which 
implies that (ImAo)^^ = has recently been given in [13]. 

Theorem 1.3 Let V G 9JI. Then Conditions Bl, B2, Al e and A2, A3 hold for all eigenvalues 
e of L$. In fact, the maps A i— > P (L\ — ie) _1 Po an d A i— » Qp.\, appearing in Conditions A2 
and A3, extend analytically (as B(H)-valued and H-valued maps) to a complex neighbourhood 
of\ = 0. 

1.2.2 Systems far from equilibrium 

We assume that condition C is satisfied. The following result is analogous to the one of Theorem 
1.2. Let P be the orthogonal projection onto the kernel of Lo, P = 1 — P, and denote by L 
the restriction of Lo to Ran P. 

Theorem 1.4 Suppose A i— » ip\ is differ entiable at A = as a map from R to H and denote 
its derivative at zero by ip' . Then we have 

lim PWP(L Q - ieytpWPipo = PWPip' . (17) 

A special case (which is of interest in concrete applications) is given by a reference state ip\ = ipo 
which does not depend on A, or by interactions satisfying PWP = 0. In either case (17) shows 
that P* G Ker A . 

2 Examples 

A detailed description of the theory of ideal quantum gases is given in [10, 25, 27, 22]. 
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2.1 Reservoirs of thermal Fermions 

The Hilbert space of states of an infinitely extended ideal Fermi gas which are normal w.r.t. 
the equilibrium state at inverse temperature < (3 < oo is 

ft R = JT_®jr_ ( 18 ) 

where T- = ©„> ^antisym^ 1 1S tne antisymmetric Fock space over the one-particle Hilbert 
space ft = L 2 (R 3 ,d 3 k) (momentum representation). The thermal creation operators (distribu- 
tions) are given by 

a}{k) = ~J^a*(k)® \ T _ +{~1) N ®^a(k), tel 3 , (19) 

whith (i,p(k) = (1 + e^^))" 1 and where 

k i — ► u(k) > (20) 

is the dispersion relation of the fermions considered. A typical example are non-relativistic 
fermions for which u)(k) = \k\ 2 . The a and a* on the r.h.s. of (19) are the ordinary fermionic 
Fock annihilation and creation operators which satisfy the canonical anti-commutation relations 

{a(k),a*(l)} = S(k-l). (21) 

The number operator N in (19) is given by N = J R3 a* (k)a(k) d 3 k. Relations (18) and (19) 
constitute the so-called Araki- Wyss representation of the canonical anti-commutation relations 
[5]. Smeared-out creation and annihilation operators are defined by a*(f) = J R3 f(k)a*(k)d 3 k 
and a(f) = J R3 f(k)a{k)d 3 k 1 for / e Sj, and where / stands for the complex conjugate. One 
shows that the fermionic creation and annihilation operators are bounded, satisfying ||a^(/)|| = 
|| /|| The von Neumann algebra 9JIr is generated by the thermal creation and annihilation 
operators {a* (/) | / G f)}. The vector 

= Q ® Q e « R , (22) 

where f2 is the (Fock) vacuum vector in J 7 ^, represents the KMS state w.r.t. the dynamics 
1 1 ^ c itLR a*(f)e- itL « = a*{e itu} f), where 

L R = dT(w) <g> Ir- - ® dT(w), (23) 

and dr(cj) = J R3 ui(k)a*(k)a(k) d 3 k is the second quantization of multiplication by oj(k). The 
action of the modular conjugation Jr on creation operators is 

Jna}(f)J R = [(-1) N ® a*{^T~^J) + a(^J) ® (-1)^ ® (-1)". (24) 

Now we describe the interaction with the small system. The latter is described at the 
beginning of Section 1.1. For m, n > 0, m + n > 1, let h m , n be maps from R 3m x R 3 ™ to 
the bounded operators on JC. For simplicity of exposition, we assume that those maps are 
continuous and have compact support. We define — (k\, . . . , k m ) € R 3m and put 

a}(k^) = a* p (k 1 )---a* p {k m ), (25) 

and similarly for ap{l^). Set 

v m , n = f dk (m) dl (n) /i roi „(fe( m \z^)®a^(fe (m) )a /5 (« ( " ) )+ adjoint. (26) 

It is well known that 

V m , n = TT(v„ hn ) e 971, (27) 
see (6). The following result follows from Theorem 1.3. 

Theorem 2.1 Suppose V is a linear combination of terms (27). Then all conditions Bl, B2, 
Al e , A 2, A3 hold. In particular, the conclusions of Theorems 1.1 and 1.2 are valid. 
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2.2 Reservoirs of thermal Bosons 

The Hilbert space describing states of an infinitely extended ideal Bose gas which are normal 
to the equilibrium state at inverse temperature < (3 < oo (without Bose-Einstein condensate) 
is 

H K =F+®F+ (28) 

where T + = ©„> ®sym¥> is the symmetric Fock space over the one-particle Hilbert space 
= L 2 (R 3 , d 3 k) (momentum representation). The thermal creation operators (distributions) 
are given by 

a}{k) = ^/l + P f } a*(k) ® % + % ® y/pp a(k), keR 3 , (29) 

where pp(k) = (eP u W — I) -1 is the momentum density distribution, and 

k i ► uj{k) > (30) 

is the dispersion relation. We shall consider for simplicity of the exposition w's such that pp{k) 
is locally integrablc in R 3 . A typical example is cu(k) — \k\ (massless relativistic Bosons). The 
a and a* on the r.h.s. of (29) are the ordinary Fock annihilation and creation operators which 
satisfy the canonical commutation relations 

[a(k),a*(l)] = 6(k-l). (31) 

The representation (28), (29) is the so-called Araki-Woods representation of the canonical com- 
mutation relations [4]. Smcared-out creation and annihilation operators are defined by a*(f) = 
J" R 3 f{k)a*(k)d 3 k and a(f) = J R3 f(k)a(k)d 3 k, for / e $j, where / stands for the complex con- 
jugate. The von Neumann algebra 9JIr is generated by the Weyl operators Wp(f) = c 1 '*' 3 ^, 
/ S f), where 

Mf) = ^%(f) + ap(f)) (32) 

is the selfadjoint field operator. The unbounded operators a^(/), ap{f) and </>/?(/) are affiliated 
with 2JIr. The vector 

<f>p = 0®Q e H R , (33) 

where is the (Fock) vacuum vector in !F + , represents the KMS state w.r.t. the dynamics 
1 1 ► e itLR Wp(f)e~ itLR = Wp{c lt "f), where 

L R = dT(w) ® - ® dl», (34) 

and dr(w) = J R3 uj(k)a*(k)a(k) d 3 k is the second quantization of multiplication by uj(k). The 
modular conjugation Jr acts on creation operators as 



Jr<4(/)Jr - ® a(^TTp^f) + a{y/ppf) ® 11^ + . (35) 

We now turn to the interaction between this reservoir and the small system (which is described 
at the beginning of Section f .1). It is given in a similar way as for Fcrmions, but, for technical 
reasons, it has to be restricted to at most quadratic expressions in the thermal creation and 
annihilation operators ("two-body interactions"). Given hi and h 2 , continuous maps from R 3 
and from R 3 x R 3 to the bounded operators on K respectively, s.t. h 2 (k, I) — h 2 (l, k)*, we set 

v= [ [h^^a^ + h^kf ®a [} {k)\d 3 k+ { h 2 {k,l)®a* l) {k)a l3 {l)d 3 kA 3 l. (36) 
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For the purpose of exposition we assume that have compact support. We introduce 

C(s) = f c 2sw W||e- sH Mfc)c sH || 2 d 3 fc 

+ / e- 2s(w(fc) - w(()) ||e-^ 2 (/M)c sH || 2 d 3 kd 3 l, (37) 

which measures the regularity of the integral kernels ft.1,2 relative to the Hamiltonian H. 
Theorem 2.2 Suppose the functions hi and h 2 are such that 

sup C(s) = C < 00. (38) 

-0/2<s<0/2 

There is a constant X (depending on (3) s.t. if\\\ < A then all the conditions Bl, B2, A1-A3 
are met. 

Remarks. 1) If dim/C < 00 then (38) is satisfied. 

2) We can replace the condition of h\$ having compact support by a suitable weaker decay 
condition. Moreover, assuming h\^ to satisfy certain regularity properties at k, I — we can 
treat discontinuous ft. 1.2. 

3) One may add to v, (36), integrals involving a l g(fc)a / g(Z), ap(k)a^(l), a^(k)a^(l). 

4) If hi = then the theorem holds for all A G R (see the remark at the end of the proof of 
Theorem 2.2, Section 4.4). 



2.3 Systems with several reservoirs 

We consider R > 2 fermionic reservoirs coupled through a small system. The Hilbert space 
is given by Ti = Hs <8> Hr x <8> • • • <£> Hr, h , with an uncoupled dynamics generated by Lq = 
Ls + £r.i + • • • + £r h - We may choose the reference state to be a product of equilibrium 
states, -00 = ® ® • • • <8> Denote by Jr^s and Ar^s the modular conjugations and 
the modular operators associated to the pairs (97t Ri , and (97ts> ^/3 S ), respectively. Take 
Vi,...,Vr to be linear combinations of the form (27), as in Theorem 2.1, representing the 
couplings to the reservoirs, and suppose that 



(A s /2 ® A^ 2 ) V, (A" 1/2 ® A Rs 1/2 ) e Wl s ® 97Tr. 



We set 



W = E - ( J s A s /2 ® ^ a r s 2 ) ^ (^sA s /2 ® J Ri A^ 2 ) } . (39) 



It is easily verified that VFV'o = and thus K\ipo = [Lq + \W)ipo = 0. Since the second 
part in the sum belongs to the commutant 971', K\ defines the same dynamics on 971 as does 
L + A(Vi + • • • + Vr). Thus condition C is satisfied and Theorem 1.4 applies. 

2.4 Klein— Gordon field for accelerated observer 

Let x = (x , x) be a point in Minkowski space-time ixl 3 (with metric signature (+, — , — , — )). 
The field operator satisfying the Klein-Gordon equation 

(□ + m 2 )ip(x) = 0, (40) 

where □ = <9 2 — A, A is the Laplacian, m > 0, is given by 

= / ^L= [e^-^a(k) + e^+^V W ] . (41) 

JR3 y/2bj(k) 
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The a # (fc) in (41) are the usual bosonic creation and annihilation operators satisfying (31), 
and 



w(fc) = (k 2 + m 2 ) 



2\V2 



(42) 



Let <S(R 4 ;R) and 5(M 4 ;C) denote the real valued and the complex valued Schwartz functions 
on R 4 , respectively. For / <G <S(R 4 ; C) we define the smeared field operators, acting on bosonic 
Fock space T+ (see after (28)), by 



<p[f} = I f(xMx)d 4 x. 
Jr 4 



(43) 



The adjoint of <p[f] is <p[f], where / stands for the complex conjugate of /. For / e 
we define the unitary Weyl operator W[f] = e lv ^. 

We introduce the "right wedge" Wr = {x e R 4 | \x°\ < x 1 }. Let 9Jl R C B(T+) be the von 
Neumann algebra generated by 



[W[f] | / G 



supp(/) C Wr}, 



(44) 



where supp(/) denotes the support of /. Elements of 9JTr are interpreted to be those observables 
which can be measured in the space-time region Wr. 

It is well known (Reeh-Schlieder) that the vacuum vector £!r e T+ is cyclic for the ^algebra 
of all polynomials in <p[f], f € 5(M 4 ; C). We may thus define an antiunitary involution Jr on 

by 

J R <p[f(x))J R = <p[f(-x°,-x 1 ,x 2 ,x% (45) 
for all / e 5(M 4 ,C). Rclativistic boosts in the x 1 -direction are given by 

B T : (x°, x 1 , x 2 , x 3 ) i— » (x° coshr + x 1 sinhr, x° sinhr + x 1 coshT, x 2 , x 3 ), 

for t G K. The action of B T lifts to Fock space according to 



ip[f o B- 1 



-WLu. 



where Lr is the sclfadjoint Liouville operator on T + given by 

L R = dr([-A + m 2 ] 1 / 4 x 1 [-A + m 2 ] 1 / 4 ). 
The map B T leaves the wedge Wr invariant so 

a T K {A) = e iTLR Ae- iTLn 
defines a group of * automorphisms of 971r. 



(46) 
(47) 
(48) 



Theorem 2.3 ([7]) The state on 9JIr determined by the vacuum vector Or e T + is a (2tt, a R )- 
KMS state. The modular operator associated to (9JIr,Or) is Ar = c - 27ri n and the modular 
conjugation is Jr. 

An observer accelerating in the x 1 -direction with constant acceleration a > 0, with position 
(1/a, x 2 , x 3 ) at the instant x° = 0, describes the curve [24] 



(49) 



" x° ' 







x 1 
x 2 




1/a 

x 2 








X 3 




X 3 



where e r{x > = ax° + ^a 2 (x ) 2 + 1. The (x') 1 -axis of the accelerated observer's instantaneous 
rest frame is given by the half line emanating from the origin and passing through the point 
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(49), in the coordinate system (x° , . . . ,x 3 ). The curve B T / X o^[0 1/6 x 2 x 3 ] T traces the point 
which is perceived by the accelerated observer to be at rest, lying on the (a;') -axis at \n(b/a)/a 
(in a natural coordinate system) [24] . 

According to (46) and (49), the field at the point of the accelerated observer at time x° is 
given by 

<p{x a ,x) = e iT< - x ^ L «(p(0,l/a,x 2 ,x 3 )e- iT< - x ^ Ll <. (50) 

Similarly, the field at the moment x° , in a space point x which is at rest as seen by the 
accelerated observer, with spatial coordinates (S',x 2 ,x 3 ) in the instantaneous rest frame, is 
given by 

<p(x°,x) = c iT ( xO)L «^(0,c s '/a,x 2 ,x 3 )e~ [T{x ° )LR . (51) 

The proper time (instantaneous time) on a trajectory of an accelerated observer in x 1 -direction 
with constant acceleration a > is r a = r/a, [24, 31]. Theorem 2.3 thus means that the field, in 
the vacuum state fijt, is perceived by the accelerated observer as a thermal field at temperature 
/3 = 27r/a ("Unruh Effect"). 

To measure this temperature concretely, the accelerated observer carries a thermometer, 
which we model by a small system, as described at the beginning of Section 1.1. The interaction 
between the thermometer and the ambient field is given by the operator [32, 11] 

V = G®1 K ®( d 4 x 5(x°) p(x)(p(x), (52) 

or by a sum of such operators, where the integral is effectively taken over the time-slice x° = 0. 
Here, G is a bounded selfadjoint operator on JC and p is a smooth function compactly supported 
in a small neighbourhood of the position of the accelerated observer at time zero, (1/a, 0, 0). 

Theorem 2.4 Suppose that sup 0<s<7r \\c~ sH Gc sH \\ < oo. Then all the conditions Bl, B2, 
A1-A3 are met for arbitrary A G R. 

We give a proof of this result in Section 4.5. Theorem 2.4 is useful in the mathematical 
analysis of the Unruh Effect, [11]. 



3 Hamiltonians with degenerate spectrum 

We consider a small system with dim/C < oo coupled to a reservoir of relativistic massless 
Bosons (dispersion relation to(k) = \k\) at temperature B > 0. The interaction is 

V = G®1 C 2 ®<i> {g), (53) 

where G is a selfadjoint operator on JC, g € L 2 (K 3 , d 3 fc), and where the thermal field operator 
<t>p{g) is defined in (32). 

Lemma 3.1 Assume that g has compact support and is continuous, except possibly at \k\ = 0, 
and suppose that 

lim MM = 7 , (54) 
IfcHo \k\P 1 ; 

for some p > — 1 and some 7 > 0. Then the system introduced in this paragraph, with interaction 
(53), satisfies all conditions Bl, B2, Al e , A2 and A3, for any value of the coupling constant A. 

We do not optimize the ultraviolet behaviour of g here. A proof of this result is easily obtained 
by repeating the proof of Theorem 2.2. 
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Theorem 3.2 Let JC = C 2 and H = diag(e, e), for some e e M. ITe Ziawe dimPo = 4. Sef 
e=||. 9 /V^H 2 /2 andr^ 2 ^. 

1. If — 1 < p < — 1/2 t/ien t/ie teweZ s/ii/i operator A does not exist. 

2. If p = —1/2 f/ien i/ie feve/ shift operator A exzsis and is given by 

Ao = (£ + «?)G 2 O 31 C 2 - (£ - ® CG 2 C - 2i?yG ® CGC, (55) 

w/iere C is the antiunitary map taking the complex conjugate of coordinates in the canonical 
basis ofC 2 . 

3. If p > —1/2 then A exists and is given by (55) with r\ = 0. 

This result is obtained by explicit calculation, see Section 4.6. Relation (55) gives 

ReA = £(G 2 <g> H C 2 - H c = O CG 2 C) (56) 
ImAo = n(G ® 1 C 2 - ]1 C 2 O CGC) 2 > 0. (57) 

For p > — 1/2 we have ImAo = 0. The particular structure of Ao allows us to find its spectrum 
easily. Let ai, ct2 be the (real) eigenvalues of G, with associated eigenvectors xii X2- Then the 
eigenvectors of A are \i ® Xj> hj G {1; 2}, and the spectrum of A is 

o-(Ao) = {0,0, z,-z}, (58) 

where z — (ai — a 2 )[^(ai + a 2 ) + irj(ai — a 2 )]. These observations prove the assertions of 
Theorem 1.1, part (c). 

A physically more interesting model is given by a small system with a non-degenerate ground 
state and two degenerate excited states. We have JC = C 3 , H = diag(e, /, /) with A = / — e > 
and dimP = 5. Take a coupling of the form (53), where 



G = 



a b 
a c 
b c 



with a,b,ce M. (59) 



Theorem 3.3 Consider the interaction (53) with G as in (59) and assume g(k) satisfies (54) 
with p > —1/2. Then A exists and has spectrum 



( IS p> JL ^ *i 

<7(Ao) - |0, 0, is(a 2 + b 2 ) coth(/3A/2), (a 2 + b 2 )(±a + - )|, (60) 

where s = ttA 2 J s2 da\g(A, a)\ 2 and a = P.V. (g, ^^ g) (principal value). 



Remark. The operator Ao is independent of c, see (108). This comes from the fact that c 
governs transitions between the two degenerate energy levels (Bohr frequency zero) and g has 
a "mild" singularity at the origin (p > —1/2). We understand this independence of the level 
shift operator of c as the origin of the degeneracy of its kernel (the coupling is "not effective" 
since one may choose c = 0). A exists and depends on c in case p — —1/2, c.f. after (109). 

The violation of the Fermi Golden Rule Condition (dimKcr Ao > 1) calls for a more detailed 
examination of the time ergodic properties of the system. For the purpose of exposition, we 
consider a translation analytic [20] model in what follows. The Feshbach map applied to the 
spectrally translated L\(t) is [8] 

F(L x (t)) = -X 2 PoIPo(L x (t))- 1 P IPo = -A 2 limP /P (L A - ie^/VP., (61) 

<40 
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where we write L instead of L . We take the convention that the spectral deformation is 
implemented in such a way that the the spectrum of L\ gets pushed into the lower complex 
plane for the range of r considered (hence the sign e J. 0) . We expand 

PoIP (Lx - ie^PoIPo = A (e) + A 2 A^(e) + 0(A 4 ), (62) 

see also (12) and (1), with a remainder 0(A 4 ) uniform in e for small e. The following result 
describes the spectrum of A + A 2 Aq, where Aq = lim e j A (e). 

Proposition 3.4 Consider the interaction matrix G, (59), with a ^ 0, 6 = 0. The operator Aq 
exists and the spectrum of A + A 2 A in a neighbourhood of order A 2 around the origin consists 
of two eigenvalues, one at the origin and another one, z, which satisfies 

Imz = a 2 c 2 ^(A)+c%, (63) 

where ^i(A),^ 2 > 0. £i(A) is of order A 2 for small A, while & does not depend on A. 

The operator A' exists also when b ^ and it is possible (albeit much longer) to calculate the 
spectrum of Ao + A 2 A in that case. We give explicit expressions for £i, £2 in equations (117), 
(118). 

It follows from (63) that the degeneracy of the eigenvalue zero of A is lifted by the pertur- 
bation A 2 Aq. By the isospectrality of the Feshbach map one concludes from (61) and (62) that 
Ker L\ is simple, provided A 7^ is small enough. Moreover, it is shown in a standard way that 
all non-zero resonances of L\ have strictly negative imaginary part. This implies the property 
of return to equilibrium with exponentially fast convergence, proportional to e~ ?A *, for some 

£>o. 

4 Proofs 

4.1 Proof of Theorem 1.1 

We realize the Hilbert space of the Gibbs state $>p of the small system by Hs = /C®/C. Let {Ei}, 
Eq < Ei < • • • , denote the spectrum of H and let ipi the normalized eigenvector associated to 
Ei. The ipi form an orthonormal basis of JC. Then 

i>0 

where Zp is a normalization factor. It is manifest that is* p = 0, c.f. (2). The modular 
conjugation Js has the explicit action 

J s iP (f) = (f) (65) 

where C is the antilinear operator whose effect is to take the complex conjugate of coordinates 
of vectors in K., relative to the basis {fi}. As is easily verified, we have the relations Js^p = * p 
and JsL$Js = —Ls (c.f. (2)). The spectrum of L$ consists of the eigenvalues {Eij :— Ei — 
Ej\i,j>0}. 

Let e be an eigenvalue of is- A basis for RanP e is given by 

{<Pi ® ifj $/3 I s.t. Ei - Ej = e}. (66) 

According to the terminology introduced in [9], Appendix B, if e is an eigenvalue of Ls, then 
{e} is called a nondegenerate set if and only if the following holds whenever i is s.t. Ei — e lies 
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in the spectrum of H: if Ei — Ej = e and Ei — Ey = e then j = f. It should be noted that {0} 
is a nondegenerate set (in the sense of [9]) if and only if all Ei are distinct, i.e., if and only 
if H has simple spectrum. If the spectrum of H is simple then a basis of Ran P is given by 
{ipi ®tpi® Qp}. 

Proof of Theorem 1.1. (a) Let pi — \ifi)(<fii\ denote the orthogonal projection onto Ctpi. We 
have 

P e = J2 P^®P3® P K■ ( 67 ) 

Ei,j=e 

Taking into account Jsfi <E> <Pj = <Pj <8>(Pi, Jr3>/3 = 3>/3 we obtain for any x^TC 

= {<Pj ®W® ®0, x) <Pj ® ifi O 4>/3. (68) 

To arrive at the second equality we use that J is antiunitary. The equivalence Eij = e ^==^ 
Ej^ — — e and (68) show that 

JP e J = P_ e . (69) 
Since JL J = — L , J 2 = 1 and J is antilinear, we have 

JP e (L e - e - ie)~ 1 P e J = -P_ e (Io " + e - ^P-e- (70) 
Assume that Al e holds. Using (69), (70) and JIJ = —I, c.f. (8), wc obtain 

JA e (e) J = -P_ e /P_ e (L ~ e + e - ie)- 1 P_ e /P_ e = -A_ e (e). (71) 

This shows assertion (a) of Theorem 1.1. 
(b) For e > we have the representation 

2Re(Io - e - ie)" 1 = i f dt e-^sgn^e-' 1 ^-^ 
Jr 



(72) 



where we set RcA = ^(A + A*) for an operator A, and sgn takes the value 1 for arguments 
t > and the value — 1 for t < 0. We thus obtain 

2ReA e (e) = i f dt e~^sgn(t) \pje-'< L ^ IP e - P e IP e IP e ] , (73) 
Jr 1 J 

where we use P e = H — P e . The second term on the r.h.s. of (73) vanishes because the integrand 
is an odd function. Using P e (L — e) = we arrive at 

2ReA e (e) = i f dt e- e|t| sgn(t)P e e itio Ie' itL ° IP e . (74) 
Jr 

The decomposition I = V — JVJ, (8), leads to four terms in the integral of (74). The contri- 
bution of the two "cross terms" where both V and JVJ occur is 

-P e e itLo JVJe- itLo VP e - P e e itL °Ve' itL °JVJP e . (75) 

The following little result is proven below. 
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Lemma 4.1 Let M be a self adjoint operator affiliated with 9JI, s.t. MP e and JMJP e are 
bounded operators. Then 

P e e itL °Me~ itL ° JMJ P e = P e JMJ e itL °Me~ itL ° P e , (76) 

for all t e M. 

Lemma 4.1 shows that e ltL oy e ~itL - m ^ e second term of (75) can be commuted through JVJ, 
so 

(75) = -P e JVJ r c -i*(^o-e) + c it(L„-e)l y ^ 

This is an even function in t so it does not contribute to the integral in (74). It follows that 

2ReA e (e) = i J & e~ e ^sgn(t)P e [ve- it(i °- e) V + JVe~ it{La+e) V j\ P e . (78) 

In order to calculate matrix elements of (78) we use the following fact. If M is any operator 
affiliated with 971 and if B € B(JC) then 1/c ® B <g> ]l^ s leaves the domain of M invariant, and 
the two operators commute strongly on that domain. (This holds since Ihc ® B ® e 971'.) 
A similar statement holds for M' affiliated with 971' and B ® 31x; <8> lw s ■ 
Let (fc,Z), (k',V) be indices s.t. Sfe,j = ,v — e. Then 

(<p k ®ipi® $/3, 2ReA e (e) <p fe , <g> y>;, <g> $g) 

= i/ di e~ e|t| sgn(i)U !r (^ fc <g> tp l <g> $0, ye _it(Lo_e V^ <g> ^ <g> 

+Sk,k' (<pv ®<Pk®&p, Ve lt{ - Lo+e ^Vipi 0^0^)1 . (79) 



Consider the case e = 0. Since all eigenvalues of 7? are simple we must have k — I and k' = I'. 
The Kronccker deltas 8ij' and 5k, k' in (79) force all indices to be the same: k = I = k' = I'. 
Thus the integrand in (79) becomes an odd function of t and the value of the integral is zero. 
This shows that ReA (e) = 0, so A (e) = iTo(e), where 

r (e) := l[Ao(e) - A (e)*] = P IP 6 P IP > (80) 
It (L y + e 2 

is manifestly a positive definite operator. (We write L for T? Q in (80).) 

Next we verify that To(e) has real matrix elements in the basis {ifi <g> tp% <8> We use 

again the properties of J, as in the proof of (a), to see that JTo(e) J = To(e) (this can also be 
viewed as a consequence of assertion (a) in the theorem and the fact that Ao(e) = iTo(e)). It 
follows that 

(ifi ® ifi ® $/3, T (e)<pj ® ipj ® $g) 



= (J tp t ® ipi ® $/5, r (e) Jfj ® ipj ® $/3) 
= ® <p» ® $ / a,r (e)^j ® <Pj ® $/?), 

so the matrix elements are real. 

We have now shown assertion (b) of Theorem 1.1, modulo the easy proof of Lemma 4.1. 

Proof of Lemma 4.1. First note that e itL °Me- itL °P e , P e e itL °Me' itL ° and P e JMJ are 
bounded. Let Xn — x(\M\ < n) be a spectral cutoff operator. Then M n := XnM G 971, and 
JM n J e 971', for all n. The l.h.s. of (76) is we weak limit of 

P e e itL °M n e- itL °JM n JP e = P e JM n Je itL °M n e- itL «P e , (81) 

as n — > oo. The equality holds since e ltio • e _I * L ° leaves 971 invariant. The r.h.s. of (81) 
converges weakly to the r.h.s. of (76), as n — > oo. ■ 
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4.2 Proof of Theorems 1.2 and 1.4 



We prove Theorem 4.2 below, which covers the proofs of both Theorems 1.2 and 1.4. 

Let L be a selfadjoint operator on a Hilbert space H and let W be an operator on H s.t. 
-ft^A = + AW defines a closed operator for X £ U, where U C 1 is a neighbourhood of the 
origin. Assume that we have K\ip\ = eip\, for some e £ R and for all X £ U, where ip\ £ Ti. Let 
P be the selfadjoint projection onto the eigenspace of L associated with e, and set P = 1 — P. 
We allow the case dimP = oo. Let L = PL P l~R an p and K\ = PK\P t Ran p denote the 
restrictions of L and K\ to RanP, respectively. We make the following hypotheses: 

HO There is a sequence £„ in the open upper complex half plane C+ converging to a ( £ C+ , 
such that all and C belong to the resolvent sets of K\ for all X £ U. 

HI PW and WP are bounded operators on H. 

H2 A i ► PWP(K\ — £„) _1 P is continuous at A = as a map from R to the bounded operators 
on H, for every ( n . 

H3 A ip\ is differentiable at A = as a map from R to H, with derivative at zero denoted 
by %■ 

Theorem 4.2 Assume hypotheses H0-H3. Then we have 

lim PWP(L - e - ie)- l PWP^ - PWP^' . (82) 

ej.0 

If HO hods for a sequence ( n £ C_ converging to a £ £ C_ then (82) holds with the limit 
replaced by lim e fo- In case HO holds for two sequences, one in the upper half plane, the other 
in the lower one, then (82) holds with the limit replaced by lim e _^o- 

Remark. One can show that under suitable regularity conditions on W, the operator 
PWP(L — e — ie)~ 1 PWP has limits as e J. and as e | 0. Theorem 4.2 does not presuppose 
the existence of these limiting operators, though. 

The proof of Theorem 1.2 is an application of Theorem 4.2, where K\ is the selfadjoint 
standard Liouvillc operator L\ (7) and where e = and tp\ = flp.\ is the interacting KMS 
state (9). HO is satisfied since L\ is selfadjoint and H1-H3 follow from Al -A3. To prove 
Theorem 1.4 we apply Theorem 4.2 with K\ given in Condition C after (11). Since W one 
easily sees that conditions H0-H2 are met, and H3 is satisfied by the assumption in Theorem 
1.4. 

Proof of Theorem 4-2. In the spirit of a Feshbach-type argument we project the equation 
(K\ — e)ip\ = onto RanP and RanP, 

PWPipx = -PWP^x, (83) 
P{K x -e)Pi> x = -XPWP^x- (84) 

Take a fixed element ( n from the sequence in HO. We add the vector (e — ( n )Pipx to both sides 
of (84) and "solve" for P>a, 

P>a = -XP(Kx - C„)" 1 PWPVa + (e - (n)P(Kx - CnY^x- (85) 
Substitution in equation (83) results in 

PWP^x - XPWP(Kx - (n^PWPiPx - (e - Q n )PWP{K x - CO^P^a- (86) 
Since Ptp = equation (86) can be written as 

- PW pt±—^ = PWP(Kx - Cnr'PWP^x - (e - ( n )PWP(K x - CO^^y^. (87) 
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The l.h.s. of (87) equals PWP^^ because we have PW^x = for all A (sec (83)). In the 
limit A — > equation (87) is 

PWP^ = PWP(L - CJ-WPVo - (e - C n )PWP(L - O^P^o- (88) 

The map z ^> PWP(L - z)- l PWPip Q - (e - z)PWP(L - z)- 1 ?^ is analytic in z e C+ and, 
by (88), it is constant on the convergent sequence Q n . We conclude that equality (88) holds 
with C, n replaced by any z G C+. Choose z = e + ie, with e > 0, 

PWP^ = PWP(L Q - e - k) _1 PWP^o + iePVKP(I - e - ie)" 1 (89) 

We now take e J. 0. Relation (82) follows since ieP(L n — e — ie) _1 P converges strongly to zero 
(irrespective of the sign of e). 

In case HO holds for a sequence belonging to C_ the same argument shows (82) with the 
limit replaced by lim e fo- ■ 

4.3 Proof of Theorem 1.3 

It is shown in [16] that Clpfi <= T>(e~P xy / 2 ) implies Condition B2. All other conditions are very 
easily seen to hold. The analytic extension of the resolvent is 

oo 

P (L X - iej-'Po = P (L - iej-'Po ]T(-A)" [l(L - ie)- 1 ]" P , (90) 

provided |A| < e/||J||. The analytic extension of the perturbed KMS state is given (modulo 
normalization) by 

fy,A = f>A)" / • • • / dft • • • d/3„ e~^ L °V ■ ■ ■ e-^ra ft0 , (91) 

where Tp >n = {(/3i, . . . , /3„) e K n | ft > 0, /?H h/3„ < /3/2}. Expansion (91) is due to Araki, 

[3], sec also [16, 10]. ■ 

4.4 Proof of Theorem 2.2 

A verification of conditions Bl, B2, A1-A3 is quite standard. We outline the main steps. One 
proves [30] (Theorem X.44) that V and JVJ are well defined symmetric operators with a core 
for essential selfadjointness given by K, ® /C <g> ® J^o, where To is the set of all finite linear 
combinations of functions in (g)™ ym C^°(R 3 , d 3 fc), with variable n ("the finite particle space over 
test functions"). An argument as given e.g. in [6] (Theorem 3.13) shows that Lq + XV and 
Lq + XV — XJVJ are essentially selfadjoint on K. £g> K. (g> To ® To- Thus, Condition Bl holds. 

Condition B2 follows from Condition A3 which we show to hold below. 

To verify Condition Al e we note that P e = XL s =e ® Pr, where XL s =e is the spectral 
projection onto the eigenspace corresponding to the eigenvalue e of is, an d Pr = 
is the orthogonal projection onto C&p. Denote by Np = N ® 11jf + + ljr + ® N, where N = 
J R3 a*(k)a(k)d 3 k is the number operator on T + . Since Np$p = the standard bound 

\\I{Np + 1)-^= C <oo (92) 

shows that Al e holds for all eigenvalues e oi Lq,. 

To verify Condition A2 we use the resolvent identity 

P IP (L X - ie^Po 

= PoIPoiLo-ie^Po + XPoIPoiLo-ier^iLx-ie^Po. (93) 
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One easily verifies that for any v > 0, 

IP(Np < v) e Ran P(Np <v + 2), (94) 

where P{Np < v) is the spectral projection of Np corresponding to the interval [0, v\. It follows 
that the norm of the second term on the r.h.s. of (93) is bounded above by C\\\e~ 2 \\P(Np < 
2)I\\ (we commute P(Np < 2) through the resolvent (L — ie) -1 )- This shows that A2 holds. 
Finally we check Condition A3. A Dyson series expansion allows us to write 

e -(S(L +XV)/2 e PL /2^ o 

= E (-?)"/ d Sl ---ds n a^(V)---a^(V)Slp fi , (95) 

n=0 K Z/ JQ<s n <-< Sl <P/2 

where flpfi is the uncoupled KMS state, (4), and where, for z € C, 

a z (V) = e [zLa Ve- [zLo . (96) 
Formula (95) holds if the series on the r.h.s. converges, [29], Appendix B.l. The relation 

e izLo a}{k)e~ izLo = e iz ^a*p(k) (97) 
together with its adjoint allows us to estimate the integrand of (95) by 

\\a^(V)a i s ^(V)---a^(V)Qp,o\\ 

= \\a^(V)P(Np<2n)^- 1 (V)P(Np<2(n-l))---a^(V)np, \\ 

< C( Sl )---C(s n )2 n n\, (98) 

where C(s) is defined in (37). By (38) we have (98) < (2C) n n\, so the series (95) converges 
provided 

C/3\X\/2 < 1. (99) 
For A satisfying the bound (99) the map A Qp,\ has an analytic extension, given by (95). 

1 /2 

Remark. If hi = then V is relatively bounded w.r.t. Np and the n\ in the r.h.s. of 
(98) can be replaced by (n!) 1 / 2 . Thus the series (95) converges for all values of A e K. ■ 

4.5 Proof of Theorem 2.4 

It is useful to pass from L 2 (R 3 ,d 3 x) to L 2 (R 3 , d^dfc^), where k 1 = (fc 2 ,fc 3 ) and fj, e R, 
according to the isometric isomorphism 



/^(2tt)- 1 / 2 I dKe-'^ K f{uj ± smhK,,k J -), 
Jr 



(100) 



where / is the Fourier transform of / and luj_ = (k + to 2 ) 1 / 2 . In the new space, the operator 
(47) becomes = dr(/i), and the interaction (52) becomes 

V = G®1k®<I>(<t), (101) 

where (j> is the usual zero temperature bosonic field operator acting on the bosonic Fock space 
over L 2 (R 3 , d^dfc^), smeared out with 



= (2-) 1/2 / 
Jr 



vfak 1 ) = (2tt) 1/2 / dKe-^p^sinh^/c- 1 ). (102) 
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One can use a standard Nelson-Commutator-Theorem argument to show that Condition Bl is 
satisfied (for all values of A). 

To check Conditions B2 and A3 we use the Dyson series (95), as for the thermal Bosons. 
The Dyson series converges due to the assumption sup 0<s<7r \\c~ sK Gc sK \\ < oo and the fact 
that e'^a G L 2 (R 3 , dpdk 1 -). The latter fact follows, via the transforms (102) and (100), from 
the Bisognano-Wichmann theorem, which asserts that <p[c/]f2 R is in the domain of the operator 
c -ttL r ^ smce — 5(x°)p(x) is supported in Wr. The convergent Dyson series defines the 
perturbed KMS vector which is entire analytic in A. In particular, Conditions B2 and A3 are 
satisfied. 

Conditions Al e and A2 are verified just as for thermal Bosons, Section 4.4. ■ 
4.6 Proof of Theorem 3.2 

Since L$ = we have P = 11c 2 ® l<c 2 ® P(i ® Pn , where is the vacuum vector in T + . By using 
the definition (12), the explicit form of V, (53), and the Araki-Woods representation (32) we 
arrive at 



Ao(e) = \p[g® % ® {a(v /1 + P 9) ® + ® a(VPff)} 
-1 C 2 ® CGC ® {ljr + ® a(^l + P 5) + ff) ® } 

/>OC 

x / di c 

Jo 



-etg-itLR 



G (g) 11 C 2 (g) {a*(Vl + P g) ® + ® a*(\/P5)} 
1 C 2 CGC <g> {11^ + ® a*(VTTp 5) + a*(VP 5) ® ^+} 



P. 



(103) 



We have represented the reslovent (Lq — ie) 1 in integral form. Taking into account e ltiR (a* (/)c 



itL E 



a* (e~ lta 7) and e" 1 *^ ( V + <g> a* (/))e 



i J °° di e et e ±1LUt = — (±w + ie) 1 , expression (103) reduces to (contractions!) 



]ljr + <g> o* (e Itw /), and the formula 



2A (e) 



G- 2 ( g 
+1 C 2 <g> CG 2 C ^.9 
-2G®CGC{g, 



1 



2ie 



w — ie + e 2 

1 



.9 

2k 



• 2 , 2 

2ie 



w 2 + e 2 



viewed as an operator on C 2 



) C 2 . Next we use the formula 
2ie 



/ dr h(r)- 
Jo r + e 



i7r/i(0), 



(104) 



(105) 



in the limit e | 0, for any bounded and continuous function h on R + . In order to calculate the 
limit e J. of the second term of the first summand in the r.h.s. of (104) we use (105) with 



h(\k\) 



W _ i 



\k\ 



2p 



da 



S 2 



|g(|fcU)| S 
|fc| 2 P 



(106) 



where use spherical coordinates in R 3 , and the factor \k\ 2 is the Jacobian. For p > —1/2 the 
function h is bounded and continuous. One proceeds similarly for the other two summands in 
(104). This shows parts 2 and 3 of Theorem 3.2. 

If — 1 < p < —1/2 then it is readily seen that J °° dr h(\k\) fc |^ £ 2 — » ioo, as e J. 0, and where 
h is given by (106). This shows assertion 1 of Theorem 3.2. ■ 
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4.7 Proofs of Theorem 3.3 and of Proposition 3.4 

Proof of Theorem 3.3. We introduce the following ordered orthonormal basis for Keris : 

{ipi <g> </?i, ip 2 ® (^2,^3 ® </>3,</?2 ® ¥>3, <^3 ® ¥>2 } , (107) 

where the ifj constitute the canonical basis of C 3 . A somewhat longish calculation (that is 
carried out as in the previous section) yields the following expression for A in the basis (107) 



A = 



2sinh(/3A/2) 



a 2 +b 2 

/3A/2 


-a 2 


-b 2 


—ab 


—ab 


-a 2 


a 2 e /3A/2 





ab /- 
2 "» 


ab~7 

2 


-b 2 







ab~7 
2 *> 


ab y- 
2 "> 


—ab 


ab f 
2 •> 


ab~7 
2 "» 


a 2 C+fc 2 C 
2 





—ab 


ab~7 
2 S 


ab f 
2 S 





a 2 C+& 2 C 
2 J 



(108) 



where s = ttA 2 / g2 dcr|<?(A, cr)| 2 and C = c^ A / 2 - f sinh(/3A/2) P.V. 
and is defined after (29). The kernel of A is spanned by the two vectors 



L0{k) = \k\ 



* = 



c ^A/2 

1 
1 






and x = 



e /3A/2 

1 - (b/a) 2 


1 
1 



(109) 



(of course, is proportional to the Gibbs state * p of the small system, (64)). 

In case g(k) satisfies (54) with p = —1/2 the matrix Ao is given by the following modification 
of (108). Add to the r.h.s. of (108) the diagonal matrix 2iTic 2 5 diag(0, 1, 1, 1, 1), where 5 = 
lim r ^ r J S 2 dcr \g(r, a)\ 2 and replace the zeroes in the matrix (108) by — 2mc 2 S 2smh (^ A / 2 ) _ 

Proof of Proposition 3.4- We write P, P instead of Po, Po in this proof. It follows from (62) 
and (61) that 



A' = lim PI P(L - ie^PIPiLo - ie)" 1 P/P(L - ie^PIP. 
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(HO) 



(For a translation analytic system this limit is easily seen to exist, compare with (61)). For 
b = expression (108) reduces considerably and one checks that the kernel of A is spanned 
by the two vectors ^ given in (109) and ty = [0 1 0]* = tps <£> ip^. We apply analytic 
perturbation theory (in A 2 ) to the matrix A + A 2 Aq. The correction of order A 2 to the doubly 
degenerate eigenvalue zero of Aq is given by the eigenvalues of the matrix 



D = \ 2 



(*,A£*) 
<*o,A(,tf) 



<*,A£,*„> 
(*o,A^* ) 



(HI) 



We have F(L\(t))Pc ^^ La+xv ^ 2 ^lp, a = for all values of A, where fi^o is the non-interacting 
KMS vector, c.f. (9). An expansion of this relation yields for the order in A 2 the equation 



AoXi 



0, 



(112) 



where 



Xl = 4 



d/3i 



d/3 2 PVe {fi 



2-/3i)L /2 



VSl0, o . 
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By taking the inner product of <J and \& with (112) we conclude that Aq*) = (\J> , A[,*) = 
0. Thus, since detTJ = and tr D = X 2 (* ,A(,* ) the spectrum of £> is {0, A 2 (W , A(,* )}- 
In order to calculate the matrix element (<3/o> Aq\I>o) we represent the resolvents in (110) by 
integrals, 

POO POO POO 

K' = -i\im dh dt 2 / dt 3 e~< tl+t ^ +t3) PIPI{s l )PI{s 2 )PI{sz)P 1 (113) 
£ lo J J J 

where 7(i) = e HLo Ie~ itL ° , si = -h, s 2 = -t\ - t 2 , s 3 = -ti - 1 2 - t 3 . Since PIP = we may 
drop the projections P in (113) except for the middle one. An explicit calculation gives 



7(s 2 )7(s 3 )#o ® $/3 



= c 



etc 



-is 2 A' 



C 



(/?3 <8> Tl - C ^ 2 <8> ¥>2 ® Tb + C</? 3 



0.0 



is 2 A' 



c 



(114) 



where Ti,T 2 , T 3 are vectors in the Hilbert space of the reservoir. Similarly, 



7(si)7* O*/3 

-isiA' 



ae 



¥>3 ® Si - C 2 (f 2 <g> <£ 2 ® 52 + C</? 3 ' 



no 



isiA' 


c 



5 3 , 



(115) 



where the Sj are obtained from the Tj by replacing s 2 by s\ and setting s 3 = 0. One checks 
that 

/•OO 

lim / dt 3 e" e ' 3 (tf ® */8, 7 7(si)P7(s 2 )7(s3)*o ® = 0, 
£ -*°7o 

so the middle P in (113) can also be dropped (we may take the epsilons multiplying t\,t 2 ,t 3 
in (113) individually to zero). The inner product ((115), (114)) equals 



a 2 c 2 (e It2A (Si, T\> + c- It2A (5 3 , T 3 » + c 4 «Si + 5 3 , 7\ + T 3 ) + (5 2 , T 2 )) . 



(116) 



A pretty lengthy calculation shows that the term in (116) which is proportional to a 2 c 2 yields 
the following contribution to (^o, Aq\I/o), 

ia 2 c 2 a(A) := (117) 

2ina 2 c 2 A 2 [ da f da' /"dr^a, r)\ 2 \g(p> , r - A)| 2 P(3 (r - A)(l + P(3 (r)), 
Js 2 Js 2 Ja 

where pp(r) = (c' 3 *' — While this contribution is purely imaginary the one coming from 

the term in (116) proportional to c 4 is not; its imaginary part is obtained again by a quite 
longish calculation and it is given by 



c% ■= 



(118) 



ttc 4 J daj da'j^ dr r 2 \g(a, r)\ 2 \g(a', r)\ 2 (2^(^(1+^(0) - 1 - 2 Pp {r)) \ 

This concludes the proof of Proposition 3.4. ■ 
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